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We review some basic features of dE / dx particle identification that are relevant to 
high energy physics tracking devices. Gas-based drift chambers as well as silicon 
trackers are discussed. 



1 Model of ionization 

One of the most successful models of ionization is the photo-absorption ioniza- 
tion modelt!l. The model assumes that the process of interest is dominated by 
single-photon exchange which allows one to relate the ionization process to the 
photo- absorption process for which the data is more readily available. Also the 
model is restricted to non-relativistic recoil electron, thus it does not include 
high-energy knock-on electrons (6 rays). In this model, the single collision 
cross section can be written as 
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where the first term results in the logarithmic rise of ionization, the second 
term is the Cerenkov term, and the last term is the Rutherford scattering off 
a quasi-free electron. The parameters arc: 
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The logarithmic rise can be seen in the first term by taking the limit of f3 — ► 1 
and low density e = t\ + ie2 — * 1: 

In 1 -> In — 3_ = l n 7 2 . ( 2 ) 

On the other hand, at large photon energy T (far above poles), the dielec- 
tric constant of material is given by the plasma energy hu> p : 



( huip ^ 2 
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where 



1^,M\) =2s.S A /|p( ff /cm3). ( 4 ) 



For a molecule, Z and A refer to those of the whole molecule; the plasma energy 
is proportional to the square root of the electron number density. Then, for 
j3 — > 1, we have a plateau at 

1 T 2 
ln — = — > ln — — (5) 

v/(l-/? 2 ei) 2 +/3 4 e| (^p) 2 
which the logarithmic rise (|2|) would cross at 

T 

7=z-- (6) 

This is for a given energy of the exchanged photon. In practice, the saturation 
energy is given by replacing the recoil energy T by the effective excitation 
energy I: 

7sat ~ r — • (7) 

TlLOp 

For elements, we have approximately i 

J(eV) - 13.5Z (Z < 14) , Z(eV) - 10Z (Z > 14) , (8) 

which is good to 20% except for /(H) = 19.2 eV and /(He) = 41.8 eV. For 
compound molecules or gas mixtures, a good approximation is the Bragg's 
additivity 

tagI=££^L, 0) 
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He 


41.8 


0.27 


154 


21/76 


Ar 


188 


0.82 


230 


32/115 


Xe 


482 


1.41 


341 


48/170 


CH 4 


41.7 


0.61 


68.4 


10/34 


C2H6 


45.4 


0.82 


55.3 


8/28 


C 3 H 8 


47.1 


0.96 


49.1 


7/24 


C4H10 


48.3 


1.14 


42.4 


6/21 



Table 1: 



where is the total number of electrons per molecule of a given element i. 
For hydrocarbons CfcH;, for example, nc — 6fc and uh = I, and gives around 
45 ± 5 eV for wide range of values for k and I. For gas mixtures, we have 
n% = Zifi where Zi is the total number of electrons per z-th molecule and 
fi is the fraction by volume. Even though the effective excitation energy in 
principle should depend on the_jtype of chemical bond, the above formula was 
found to work for most cases, u The heavier the element, and the lower the 
density, the higher is the saturation point 7 sat . Table [l] gives the saturation 
point for some gases as well as the saturation momentum for 7r (P^t) and K 
(i^rt). For a given type of gas, the Kit separation starts to degrade as the 
momentum goes higher than P^ at and becomes useless around P s ^ t . 

2 Mean dE/dx vs (3j 

When averaged over the recoil energy, (Q) eventually leads to the (truncation- 
improved) Bethe-Bloch formula: U 

dE 1/1 2m e /3 2 7 2 r /? 2 / T \ S\ 
^° C ^U IOg— P t( 1+ 7^)-2J ^ 

where To = rnin(T cut , T max ), T cut is the effective cutoff for the knock-on elec- 
tron energy, and T max is its kinematic limit: 

2( 7 /3) 2 m e 

where x — m e /M with M being the mass of the projectile. We have T max ~ P 
for 7a; ^> 1, or the recoil energy can be as large as the energy of projectile itself. 
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Such large recoil electrons should clearly not be included in the evaluation of 
dE/dx measured along the projectile since they will form separate tracks. The 
truncation T cut in the expression above is designed to take care of such effect. 
Typical value for T cut is 10 to 100 keV depending on the selection criteria used 
for drift chamber hits. A knock-on electron of 1 MeV will curl up with r = 1.5 
mm in a 2 Tesla field, and will probably be rejected by tracking code. Also, 
usually a fixed fraction of high-side (or sometimes also low-side) tail is rejected 
in forming the dE/dx of a track (the truncated mean method) which also leads 
to an effective T cut . Figure plots Q for argon with no T cut , T cut = 100 keV, 
and = 10 keV. The absolute value of the total relativistic rise is reduced 
by about factor of two upon applying T cut , but is relatively insensitive to the 
values of T cut in the range of 10-100 keV. 

3 Resolution 

Empirically, the resolution of gas-sampling device can be approximated by@ 

- = 0.41n-°- 43 (xP)-°- 32 , (12) 
A* 
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Belle 


52 


1.5 


1 atm 


6.6% 


5.1% (/x) 


CLE02 


51 


1.4 


1 atm 


6.4% 


5.7% (ji) 


Aleph 


344 


0.36 


1 atm 


4.6% 


4.5% (e) 


TPC/PEP 


180 


0.5 


8.5 atm 


2.8% 


2.5% 


OPAL 


159 


0.5 


4 atm 


3.0% 


3.1% (/z) 


MKII/SLC 


72 


0.833 


1 atm 


6.9% 


7.0% (e) 



Table 2: Measured resolutions are compared to the expected dE/dx resolutions that are 
estimated by the empirical formula given in the text. 



where n is the number of samples, x the sampling thickness in cm, and P 
the pressure in atm. This is obtained by a fit to actual resolutions measured, 
and includes optimizations of resolution using such techniques as the truncated 
mean method. Another method is to employ a maximum likelihood fit to all 
the pulse-heights along a track. Such technique is known to result in similar, 
if slightly better, resolutions. Table || shows a comparison of some resolutions 
measured by recent experiments and those 'expected' from the formula above. 
The formula works well except for CLE02 and Belle for which the gas contains 
a large fraction (50%) of hydrocarbon (ethane) Hydrocarbons typically have 
better resolutions for a given sampling condition n, x and P since they have 
relatively small values of mean ionization energy with respect to the number 
of electrons per molecule resulting in a larger total number of ionization. 

If each sampling is independent and there are no other sources of error 
such as electrical noises, then the resolution would scale as n~ 5 . Since the 
power on n is larger than that on x, for a fixed total length one obtains a 
better resolution for a finer sampling. At some point, however, the number of 
primary ionization which is approximately given by[] 

rip ~ 1.5Z/cm (Z : whole molecule) (13) 

becomes of order unity and there will be no gain in resolution thereafter. 
Typically, the critical sampling size is about a few mm. 

We also notice in Table || that higher pressure is effective in improving 
the resolution. As we have seen in (^|) and (Eh, however, the saturation energy 

a A better fit can be obtained by replacing (xP) with 6.83n e a;(cm)P(atm)//(eV) where n e 
is the number of electron per molecule. See Reference 1. 

'"This formula is usually good to 20%. Some exceptions are 5.2 for H2, 5.9 for He, and 44 
for Xe. 
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Figure 2: K/tt separation (number of sigma) for argon as functions of logP(GeV/c). Two 
curves correspond to 1 atm and 5 atm. 

is proportional to l/*Jp or equivalently to 1/yP and thus dE/dx saturates 
at lower energies for higher pressures. Also, a high-pressure chamber adds 
problems to design and fabrication of the tracking device, and also adds more 
materials compared to non-pressurized chambers. To study particle type sep- 
aration, we take two reference points for argon: 4.5% resolution at 1 atm and 
2.5% at 5 atm. These are quite realistic values. The resulting K/tt sepa- 
ration defined as {dE/dx\ n — dE/dx\K)/o'dE/dx is shown in Figure |] for the 
two cases. For 1 atm, 2a separation can be achieved for p < 0.8 GeV/c and 
1.75 < p < 65 GeV/c. On the other hand for 5 atm, 2a separation can be 
achieved for p < 0.9 GeV/c and 1.25 < p < 50 GeV/c, and Aa separation can 
be achieved for p < 0.6 GeV/c and 1.75 < p < 30 GeV/c. 

4 Silicon trackers 

Silicon has p = 2.33 gr/cm 3 and / = 173 eV. As a result, the saturation occurs 
at a low value of 7 sat = 5.6 which can be compared to the minimum dE/dx 
location of 7 ~ 4. There is essentially no logarithmic rise which is a general 
feature of solid or liquid, and thus K/tt separation in the logarithmic rise region 
is not practical. It can, however, be used below the minimum ionization region. 
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Achievable resolution is limited by the number of sampling even though 
the total ionization energy in 5 layers of 0.3 mm thick silicon is about the same 
as that in 1.5 m of argon gas. Since each sampling is 'thick', _the fluctuation 
in the Landau tail degrades the resolution. A Babar study □ based on the 
Vavilov model of ionization showed that 10.4 % resolution can be obtained 
by 5 layers of silicon strip detectors (0.3 mm thick) for pions at 450 MeV/c. 
Also an ALICE study based on 2 layers of silicon strip detector and 2 layers of 
silicon drift detector resulted in 10.6 % resolution for kaons at 0.98 GeV/c. In 
both cases, the numbers correspond to the truncated mean method where two 
highest pulse heights were discarded. Such resolution can give 2tr separation 
for p < 0.65 GeV/c, and quite adequate for search of new heavy charged 
particles such as stau in the gauge-mediated SUSY breaking models. It should 
be noted, however, that in order to fully utilize the low energy region where 
dE/dx rapidly rise as l//? 2 , the dynamic range of the readout system should 
cover up to around 20 times the minimum value of dE/dx. 

5 Summary 

The logarithmic rise of dE/dx and the saturation energy of the rise are well 
understood; in particular, the saturation point is higher for heavier element and 
pressurization results in a lower saturation. The measured dE/dx resolutions 
can be well approximated by an empirical formula which can be used as a 
guide to select the type of gas. The resolution is usually better for gases 
with higher content of hydrocarbons, but they typically have lower saturation 
energies. The larger is the number of sampling, the better is the resolution, 
but up to the sampling size of a few mm. For silicon trackers, the logarithmic 
rise is essentially absent, and thus only useful region is below the minimum 
ionizing point. Still, the achievable resolution is about 10% which is adequate 
for detection of new heavy particles such as stau in the gauge-mediated SUSY 
breaking models. 
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